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ABSTRACT 

Using the latest physical modeling and constrained by the most recent data, we develop a phe- 
nomenological parameterized model of the contributions to intensity and polarization maps at mil- 
limeter wavelengths from external galaxies and Sunyaev-Zeldovich effects. We find such modeling to 
be necessary for estimation of cosmological parameters from Planck data. For example, ignoring the 
clustering of the infrared background would result in a bias in of 7a. We show that the simultaneous 
marginalization over a full foreground model can eliminate such biases, while increasing the statistical 
uncertainty in cosmological parameters by less than 20%. The small increases in uncertainty can be 
significantly reduced with the inclusion of higher-resolution ground-based data. 

The multi-frequency analysis we employ involves modeling 46 total power spectra and marginaliza- 
tion over 17 foreground parameters. We show that we can also reduce the data to a best estimate of 
the CMB power spectra, and just two principal components (with constrained amplitudes) describing 
residual foreground contamination. 

Subject headings: cosmological parameters — distance scale — large-scale structure of universe 



1. INTRODUCTION 

The cosmic microwave background is arguably the 
most powerful probe of the parameters and viability of 
cosmological models. With temperature measurements 
on angular scales larger than a third of a degree ah eady 
at the cosmic variance limit (]Komatsu et al. 2010), fur- 
ther progress now depends on improvements at smaller 
angular scales as well as improved measurements of the 
polarization at all angular scales. Advances in tempera- 
ture anisotropy measurements come most recent ly from 
the At acama Cosmology Telescope (ACT; Dunkley et al. 
(2010)), soon to be followed by new results from the 
Sov^ Pole Telescope (SPT; Keisler et al. 2011). By 
the end of 2012 we expect a dramatic improvement from 
Planck, which is likely to leave very little room for further 
improvement in the measurement of the primary CMB 
temperature power spectrum. These improved measure- 
ments will translate into much tighter constraints on cos- 
mological models (e.g. Pl anck Collaboration||2006| . 

At angular scales smaller than a tenth of a degree, 
extragalactic foreground^ become important for three 
reasons: 1) the CMB power spectrum is dropping in am- 
plitude, 2) cosmic variance is smaller and 3) foregrounds 
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are growing in amplitude. At sufficiently small angu- 
lar scales, foregrounds become the dominant signal at 
all CMB frequencies. Furthermore, unlike galactic fore- 
grounds, they are statistically isotropic and thus cannot 
be avoided by masking regions of higher contamination. 
Their modeling is an unavoidable necessity. 

In this paper we present a parameterized, physically- 
motivated, phenomenological model for the extragalactic 
foregrounds and consider it in the context of extracting 
cosmological parameters from the primary CMB anisot- 
ropy. We demonstrate that for an analysis of Planck 
data, such modeling is necessary to avoid significant 
biases in cosmological parameter estimates, but that 
marginalization over even a very rich foreground model 
is essentially "for free"; the foregrounds are sufficiently 
orthogonal to the primary CMB that the statical errors 
on cosmological parameters are degraded by at most 20% 
for Us and less than 10% for other parameters. With the 
addition of higher resolution ground-based data or non- 
CMB Planck bands to clean the foregrounds, the degra- 
dation is reduced to a few percent for all parameters. 

The importance of extragalactic foregrounds for CMB 
analysis has been recognized for a loiig tirne ([Tegmark 



fc Efstathioull995!| Bouchet fc Gispert|1999||Kjiox|ltj001 

Tegmark ct al. 20goP ' ' 
2008 Dunkley et al 



ILeach et al.||2008[ [Cardoso et al.| 



2010). Potential biases from extra- 



galac tic contan iinams have been pointed g u t previously 



by Kn ox et al. (1998); Santo s et al. 



( |2005| ; |Serraet al.i2008] ) and |Taburet 
tmguishmg our work is the siniultaneoi 



(|2003[): 
it et al. 



Zahn et al. 
(2009|). Dis 



aneous consideration of 
all foreground components necessary for an analysis of 
Planck data, and physical modeling of these components 
informed from r ecent m easurem ents beyond t he damp- 
ing tail by SPT |Hall et aL,,2010| jVieira et al.i2010, ,Shi^ 



rokoff et al.|2010[ ) and AC'l ( [Dunkley et ai.|20in| rm tins 
paper, we will consider the foreground power contribu- 
tions from shot noise due to radio galaxies and dusty star 
forming galaxies (DSFGs), the clustering of the DSFGs, 
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the thermal and kinetic Sunyaev-Zeldovich effects (tSZ 
and kSZ), and correlation between the tSZ and DSFG 
components. We now turn to summarizing recent devel- 
opments in both modeling and measurements of these 
extragalactic foregrounds. 

Our understanding of the power spectrum due to DS- 
FGs at frequencies relevant for CMB analysis has been 
rapidly improving. We demonstrate here that for analy- 
sis of Planck data, the effects of DSFG clustering are the 
most important of the foregrounds to model. Although 
it is the most important effect, it has been almost en- 
tirely ignored by previous cosmological parameter error 
forecasting work. To date, the only papers to consider 
the impact of DS FG clustering on cosm ologi cal parame 



ter estimates are |Dunkley et al. (20101 and Serra et al. 
([2008). 

DSFG clustering powe r was first detec ted at CMB 
frequencies by the SPT (Hall et al. 2010), with subse- 
quent confirmation an d irnproved constramts from ACT 
(iDunkley et al.||2010| and SPT dShirokoff et~ aD'20T0). 
I'he recent su ite of early Planck papers ( Planck Collabo- 
ration|201l"r in particular) have also provided significant 
constramts on both the amplitude and shape of the clus- 
tering power. The Planck measurements rule out many 
otherwise viable models which generally predict higher 
power (on the scales relevant for analysis of the primary 
CMB power spectrum) than observed. 

Radio galaxy source counts from high-resolution 
ground-based data are particularly useful for Planck 
since they are sensitive to the decade in brightness below 
Planck^s flux cut. The radio sources in this brightness 
range create the dominant source of shot noise power in 
most of the Planck frequencies which contain significant 
CMB informa tion. SPT meas urements of point source 
populations ( Vieira et al.|2010 ) have offered valuable in- 
formation about the amplitude of Poisson power, as well 
as the coherence of these shot-noise fluctuations from fre- 
quency to frequency. 

Recent data, as well as recent theoretical develop- 
ments, inform our modeling of the power spectrum of the 
tSZ effect — a spectral distortion that arises due to inverse 
Compton scattering of CMB photons off the hot electrons 
in groups and clusters. The magnitude of the tSZ sig- 
nal is proportional to the thermal pressure of the intra- 
cluster medium (ICM) integrated along the line of sight. 
Upper limits on the amplitude of the tSZ power ( set b; 



Lueker et al. (2010), confl rmed by 



and further tightened by Shirokoit et al.l |2U10p ) were 



Dunkley et al.|(2010| 



found to be surprisingly low compared to p redictions 
from halo model calculations (Komatsu & Seljak 2002) 
and non-radiative hydrodynamical simulations ( |White| 



et al. j200 2). Recent work has demonstrated that the 
inclusion of a significant non-thermal contribution to the 
total gas pressure in groups and clusters in analytic mod- 
els can significantly r educe the predicted amplitude of th e 
tSZ power spectrum ( |Shaw et al.|2010 l |Trac et al.|2010[ ). 
Non-thermal pressure, sourced by bulk gas motions and 
turbulence, reduces the thermal pressure required to sup- 
port the ICM against gravitational colla pse and thus the 
amplit ude of the tSZ signal. Similarly, |Battaglia et al.| 
(2010) demonstrated that the inclusion of radiative cool- 
ing, star formation and AGN feedback in hydrodynami- 
cal simulations substantially lowers the tSZ power com- 
pared to simulations that omit these processes. Current 



predictions for tSZ power from models and simulations 
are consistent with the upper limits derived from obser- 
vations. 

These recent modeling developments are supported by 
data from Planck; when the models are used to extrap- 
olate from X-ray measurements to a predicted tSZ sig- 
nal, the predictions agree with Planck SZ observations. 
Agreeme nt is seen both in observation s of single galaxy 
clusters ( Planck Collaboration|[2011h|i I and via a stack- 



ing analaysis over a broad range m X-ray lui ninositv 



down to masses as sm all as M500 
Collaboratio7I|[20Tlj| ). 



5 X IO^^Mq (Planck 



Current data provide no direct lower limits to the am- 
plitude of tSZ power due to a degeneracy w ith the kinetic 



SZ power spectrum (Lueker et al. 20101. The kinetic 



SZ effect arises due to the Doppier 'I'homson scatter 
ing of CMB photons off of regions of ionized gas with 
bu lk peculiar velocities . Upper limits on kSZ power set 
by [Lueker et al. (|2010|) an d now substantially tightened 
by ]Shirokott et al.| ( |2()To| are ruling out some models 
of patchy reionizatiom It is useful to decompose the 
kSZ power into contributions arising from an inhomoge- 
neous transition from a neutral to ionized inter-galactic 
medium, so called "patchy reionization," and those from 
the post-reionization era, the "Ostriker-Vishniac" (OV) 
effect. The former is much more uncertain than the lat- 
ter, and our best knowledge of its amplitude comes di - 
rectly from the upper limits in [Shirokoff et ah] (2010). 
The OV power level has a currentTheoretical uncertainty 
that we estimate to be about a factor of 2. Despite its 
low levels, kSZ power is a worrisome source of poten- 
tial bias of cosmological parameters since its spectral de- 
pendence is the same as the primary CMB temperature 
anisotropics. 

We expect that the only potentially significant ex- 
tragalactic contributions to polarization anisotropy are 
Poisson power from radio sources and DSFGs. A po- 
larization analog for DSFG clustering could only arise 
due to (unexpected) correlations between galaxies in the 
polarization orientations of their emission. Polarization 
signals arise from scattering off of electrons in clusters 
and groups (ISazonov fc Su n^^^aev 1999| Carlstrom et al. 



2002 



Amb lard & White 20(15 ) and i n reionized patches 
( Knox et^.||1998j|Santos et al.|2003 l, but these are also 
expected to be negligibly small" 

In addition to developing and exploring the implica- 
tions of an extragalactic foreground model that takes 
into account recent developments, we introduce a new ap- 
proach to analyzing the multi-frequency data. We show 
how the complexities of our modeling can be reduced to 
a fairly simple description of the contamination of the es- 
timates of CMB power spectra. The contamination can 
be described by just a few principal components whose 
amplitudes are constrained by CMB-free linear combi- 
nations of the auto and cross-frequency power spectra. 

The outline of our paper is as follows. In Sec. [2] we de- 
scribe our foreground models before describing our fidu- 
cial models and surveys in Sec. [3] and [4] respectively. In 
Sec. [5] we describe our general methodology before de- 
tailing our principal component approach in Sec. [6j We 
finally present our results in Sec. [7] and discuss them in 
Seel 
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2. MODELING 

2.1. Emission from External Galaxies 

External galaxies in the bands of interest are well ap- 
proximated by power-law intensities 1,^ (x v", and divide 
fairly cleanly into those with spectral indice s a < 1 (ra- 
dio ga laxies) and those with a > 1 (DSFGs) ( Vieira et al. 
20T0] ). 

External galaxies lead to anisotropy via their discrete- 
ness, usually modeled with a Poisson distribution, and 
also via correlations due to their tracing of the large- 
scale structure. The Poisson fluctuations are important 
for both radio galaxies and DSFGs, while clustering i s 
only significant for the dusty galaxies ( Hall et ar]|2010 ). 

The Poisson contribution depends on the brightness 
function, dN/dS, via 



(1) 



where Sc is the flux cut; map pixels with sources with 
S > Sc are masked. Clustering power, in contrast, scales 
approximately with the square of the mean intensity, , 
with 

Although radio sources do cluster, their mean intensity 
at the relevant frequencies is much smaller than for the 
DSFGs; sufficiently smaller that their clustering power is 
negligible. 

2.1.1. Radio Galaxies 



From Vieira et al. 



( 2010 ) we know the radio galaxies at 
150 GHz and 22(JGHz, at flux de nsities below 100 mJy , 
are described quite well by the de Zotti et al. (20051 
modej^ This model has a brightness function that is 
approximately a power-law SdN/dS oc S^'^ . This trans- 
lates into Poisson power which depends on the flux cut 
via Ce oc SJ^"*"^. Due to the inhomogeneity of the Planck 
sky coverage, Sc will vary significantly across the sky. So 
that these angular variations can be taken into account, 
we chose to model the radio galaxies in terms of dN/dS 
rather than C^. 

For frequency dependence, we assume a power law 
spectrum for each source with spectral index a = a + da 
and the departures from the mean (Sa^) = uncorre- 
lated from source to source. With these assumptions our 
power spectra from radio sources are given b}|^ 



Cf = C 



R.O 



vv 



(3) 



2.1.2. Dusty Star-Forming Galaxies 
For the DSFGs, even the contribution to the shot noise 



the Poisson power is nearly independent of fiux cut and 
we choose to build our model in rather than dN/dS. 
In that case, the DSFG Poisson contribution is given 
simply by 



C, 



D 



c 



D,0 



(4) 



A number of authors have considere d the clustering of 
the in frared background, starting with^Bond et al. I (119 861 
Further theoreti cal investigation (Scott fc While] 



19911 



Haiman fc Knox| |2000 ) was stimulated by t he de 



tection ot the mtrared ba ckground in COBE data (Puget 



et al. 1996 Fixsen et al.|1998), and the dete ction of bright 



"sub^ millimeter" galaxies in SCUBA data (Hughes et al 

19981. Subsequ ently, the clustering has been detected 

at 160 microns fLagache et al. 2007), at 250, 350 and 

500 microns by the B alloon-borne Larg e Aperture Sub- 

miUimete r Telescope ( fViero et al.||2009[ BLAST) and at 
^-^ — 



217GHz (Hall et al. 



2010| [Dunkley et al.||20l0t. Re 

Cc 



cent Planck measurements of the CIB ( jPlanck CoUab 
orationi |2011r|) have extended to much larger angular 



scales than before at 217 GHz, 353 GHz, 545 GHz, and 



857 GHz a nd recent Herschel measurements (Amblard 
et al.|20lT I have tightened up the BLAST measurements 
and extended them to smaller angular scales. The field 
is rapidly evolving 



For the cl ustering, we assume the same model as in |Hall 



et al. (20101, extended to phenomenologically include the 



consequences of non-linear clustering by including a mul- 
tiplicative factor which is a power- law in £ for i > 1500. 



Additionally, the Hall et al. (2010) model has a spectral 
index that varies with because t 
tion varies with £. 



ihe radial window func- 
This effect leads to corrections on the 
order of only 1% across the relevant frequency range so 
we ignore it. Thus, the DSFG clustering power spectra 
are given by 



2 



ac 



1500 



nc 



I < 1500 
e > 1500 



falls well below the flux threshold ( [Hall et al.|2010| , thus the model by paiman ^rknox\ ( |2000[ ) 

Though they nave similar shapes in tn 



^ \Ve also k now from recent Planck results |Planck Collab<>| 
|ration|[2011m[ | that at brighter flux densities the deZiOtti model 
sigmhcantly over predicts the number counts. 

In deriving this form we have used the identity that for a zero- 
mean Gaussian random variable x, (exp(— x^)) = exp(— (x'^}/2). 
This identity and its applicability in this context, was pointed out 
to us by Challinor, Gratton and Migliaccio. 



(5) 

where <i>^^° is the Hall et al. (2010) clustering template. 

Though the same sources generate both the Poisson 
power and the clustering power, they are weighted dif- 
ferently, thus for our baseline model we conservatively 
assume no relationship between the clustering spectral 
index and the Poisson spectral index. 

To gain some idea of the range of possible shapes of the 
DSFG clustering power spectrum, we show a sampling of 
power spectra from models in the literature in Figure [ij 
They are all normalized at ^ = 3000 to highlight simi- 
larities /diff erences in shape. T he models are the fiducial 
model from|Righ i e t al. J 2008 ), the (3 = 0.6 model from 
Amblard fc Cooray' m m~and a non-linear version of 
rfnoxl (|2000|, hereafter HKOO. 

le linear regime at 



Though they have similar shapes m t 
large scales, then turn to a power-law behavior at high 
i, they resu l t from very different modeling assumptions. 
|Righi et al. ( 2008 ) associate the sour ces of infrared light 
with starbursts triggered by mergers. [Amblard fc Cooray] 
( [2007|) incorporate nonlinearities using a halo model. For 
theTlKOOnonlin' curve, we used the luminosity densities 
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Fig. 1. — Three model DSFG clustering auto-spectra at 217 GHz 
(black), and approximations to them with our parameterized model 
(solid, red), all normalized (with one exception) at £ = 3000. Our 
fiducial model is the thickest curve. Also plotted are estimates 
of the clustern g power from Planck ( Planck Collaboration. 2011r) 
and SPT (Hall et al.|2010| |. For both sets ot data points we have 
subtracted estimates ot tiie Poisson power from the reported total 
GIB power. The lowest amplitude solid (red) curve is the result of 
a "by-hand" adjustment of our model parameters to fit the Planck 
and SPT data. 

for the fiducial model of HKOO, assumed light is a biased 
tracer of mass, and calculated the no n-linear mass power 
spectrum using the prescription by 'Peacoc k fc Dodds] 
pi)96j . We will show that our phenomenological model 
is sufficient to reproduce these shapes with enough accu- 
racy for Planck cosmological parameter estimation. 

One result of the Planck measurements, available only 
after our calculations for this paper were completed, is 
that the CIB power spectrum uncertainty at £ < 2000 is 
now much smaller than before. At least two of the three 
models shown in Figure [T] that guided our understanding 
of the range of possible amplitudes have shapes that are 
inconsistent with the combined Planck and SPT data. 
That range of possible amplitudes is now given by the 
Planck CIB power spectrum measurement uncertainty. 

2.1.3. Polarization 

We expect polarized emission from the sources we con- 
sider to be very small and uncorrclatcd from source to 
source. For a collection of sources with polarization frac- 
tion /, contributing a Poisson temperature power spec- 
trum of Cn we have 



BB _ 



^TT,P 



TE _ 



TT.P 



(6) 



We parameterize both radio source and DSFG contribu- 
tions with the above forms, with / = f]j for DSFGs and 
f = fn for radio sources. 



2.2. Thermal SZ Effect 

The thermal SZ effect is a distortion of the CMB 
caused by inverse Compton scattering of CMB pho- 
tons off electrons in the high temperature plasma within 
galaxy clusters. To first order, the temperature change 
of the CMB at frequency v is given by i^T /TcMsixy) = 
f{x^)y, where /(a;^) = a::i.(coth(a;^/2) - 4), = 
hiy/ksTcMB, and y is the dimensionless Compton-y pa- 
rameter 

where the integral is along the line of sight. Tcmb is the 
CMB temperature, rig and Tf, are the number density 
and electron temperature of the ICM, respectively. 

The thermal SZ power spectrum can be calculated by 
simply summing up the squared, Fourier-space SZ pro- 
files, y, of all clusters: 



tsz 



dz- 



dV 

dz 



dlnM 



(8) 

where V(z) is the comoving volume per steradian and 
n{M, z) is the number density of objects of mass M at 
redshift z. For the latter we use the fitting function of 
Tinker et al. (2008). y{M,z,r) is the projected radial 
SZ prohle for a cluster of mass M and redshift z. Note 
that this ca lculation assumes that halos a re not spatially 
correlated; Komatsu & Kitayama ([1999') demonstrated 

usFe 



that for £ > 1000 the two-halo (or clustered) contribution 
to the tSZ power spectrum is several orders of magnitude 
smaller than the Poisson contribution given by Eq. [Sj 

To calculate the thermal SZ signal we a dopt the an- 
alytic intra-cluster gas model presented in |Shaw et al.| 
(20101. This model provides a prescription for calcu- 
lating the compton-y (or equivalently, thermal pressure) 
profiles of hot gas in groups and clusters. The model 
assumes that gas resides in hydrostatic equilibrium in 
the potential well of dark matter halos with a polytropic 
equation of state. The dark matter potential is modeled 



by a Navarro- Frenk- White profile (Navarro et al 
concentration relation 



19971 



using the h alo mass - concentration relation ot 
et al.| (2008). The model includes parameters to account 
for gas heating via energy feedback (from AGN or super- 
novae) plus dynamical heating via mergers. The stellar 
component of the baryon fraction in groups/clusters is 
determined using the stellar mass fraction - total mass 
relation observed by (Giodini et al. 2009). A radially- 



dependent non-thermal pressure component of the gas is 
incorporated by calibrating off the non-thermal pre ssure 



profile s measured in hydrodynamical simulations (Lau 
et al.| [2009). In total the model has four free parameters 
relating to astrophysical processes in groups and clusters. 
Shaw et al., (2010, ) explored the range in which these pa- 
rameters reproduce radial profiles and scaling relations 
derived from X-ray observations of nearby groups and 
clusters. 

To allow the astrophysical uncertainty to be marginal- 
ized over quickly in our MCMC chains, we perform a 
principal component analysis (PGA) described in Ap- 
pendix[B| A suite of 10,000 simulated power spectra were 
createoTeach time randomly sampling from the input as- 
trophysical parameter distribution (with the cosmologi- 
cal parameters fixed to their fiducial values described in 
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Fig. 2. — Comparison between recent models and simulations 
of the tSZ effect (black lines) and fits of our PCA model to each 
(thin red lines). The thickest red line shows the fiducial tSZ power 
spectrum used in this work. All results are plotted at f46 GHz 
and arc scaled to erg = 0.8. The blue arrow shows the SPT 95% 
con fidence upper limit on thermal SZ power at £ = 3000 | |Shirokoff| 
[eral.|2010| ). 



Fig. 3. — Comparison between recent simulations of the kSZ ef- 
fect (black lines) and our model. We consider contributions from 
the post-reionization kSZ effect (solid red line) and from patchy 
reionization (dashed red line). Note that the simulations plotted 
here assume homogeneous reionzation and thus do not include a 
patchy contribution. The SPT 95% confidence upper limit for the 
kSZ power at £ = 3000 is 6.5 fiK'^. 



Sec. [3]). We find that two principal components are suffi- 
cient to achieve 1% accuracy on the model power spectra. 

In Figure[2]we plot the thermal SZ power spectrum pre- 
dicted by a number of recent simulations (black lines) as 
well as a fit to each with our PCA model (red lines). The 
dotted line represents the thermal SZ power spectrum 
measured from the Mare-Nostrum simulation - a non- 
radiative simulation run using the smoothed-particle hy- 
drodynamics code, Gadget-2. The black solid line shows 
the results o f the non-radiative simulation of |Battaglia| 
et al. (2010) and the black dashed line the results of a 
rerun of this simulation including radiative cooling, star- 
formation and energy feedback. The dot-dashed line 
show s the 'standar d' tSZ model from the simulations 
o ^Trac et al.| (2010). The thickest red lines represents 
our fiducial thermal SZ model in this work. The blue 
point with errorbars show the recent SPT constraint on 
the amplitude of thermal SZ power a.t £ — 3000. All 
models are plotted at 146 Ghz and have been scaled to 
our fiducial cosmology. 

Our PCA model can accurately reproduce all the sim- 
ul ations in Fig [2] other th an the non-radiative simulation 
of Battaglia et al. (20101 which peaks at much smaller 
ang ular scales than the other simulations. We note that 
the Shaw et al.] (|2010|) model inherently assumes that 
some fraction of cluster gas has been converted to stars, 
whereas this simulation did not include these processes. 
Turning off star formation produces a p owe r spectrum 
that peaks at smaller scales. In Section |7.2| we investi- 
gat e the bias on measured cosmological parameters when 
the lBattagha et al.| ( |20T o') non-radiative template is used 
for the tSZ signal. We find that the PCA will adapt suf- 
ficiently to prevent a bias in the measured cosmological 
parameters. 

The final step is to determine the cosmological scaling 
of the power spectrum of our fiducial model so that the 
amplitude can be scaled accordingly in our analysis. We 



find that the tSZ power spectrum is principally sensitive 
to flm, and (78, with a particularly strong depen- 

dence on the latter. To determine the scaling we simply 
evaluate our model varying each cosmological parameter 
in the range ±25% of its fiducial value while holding the 
other three fixed (at their fiducial value). We then fit to 
the resulting power spectra, with our results summarized 
in Table □ 

2.3. Kinetic SZ Effect 

The kinetic SZ effect is a temperature anisotropy that 
arises from the Compton scattering of CMB photons off 
of electrons that have been given a line-of-sight peculiar 
velocity by density inhomogeneities in the matter field. 
We break up the kSZ into contributions from the post- 
reionization period and from a period of inhomogeneous 
"patchy" reionization. 

2.3.1. Ostriker-Vishniac Effect 

When the density fluctuations which source electron 
velocities are in the linear regime the effect is kn own 
as the Ostriker-V ishniac (O V) effect, as der ived in Os- 
triker fc Vishniac| ( |T986| ) and |Vishniac| ( |1987[ ). The post- 
reionization kSZ effect can then be modeled as the non- 
linear extension of the OV effect as described below. 



We follow the analytic prescription given in |Hu ( 2000 1 
which describes the angular power spectrum of the linear 
Vishniac effect as 



Cf = 



2^5 



dxD\ I g 




(9) 



where x is the conformal time, G is the cosmological 
growth function, Da is the comoving angular diameter 
distance, g is the visibility function, is the linear the- 
ory baryon density power spectrum and ly represents the 
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mode coupling of the linear density and velocity fields: 



with 



(10) 



/i = k • ki 

yi^ki/k 

2/2 = fc2/fc = 



1 - 2fj,yi 



yl 



(u) 



Due to an incomplete treatment of the effects of pres- 
sure feedback from baryons we slightly ov er-p r edict the 
power on very small scales. As described in Hu ( 2000 ), in 
this formulation we can consider the kSZ effect to be the 
nonlinear extension of the linear Vishniac effect. This ap- 
proximation requires replacing the linear density power 
spectrum in Eq. l9] with its nonlinear extension while 
leaving the contribution from the velocity power spec- 
trum unchanged; 



c: 



kSZ,OV _ 



2£5 



dxD\ 9 



Oh Ob 



(12) 



For the no nlinear power spe ctra (NL) we utilize the 
HALOFIT ( [Smith et al.|2003[ ) model. In this calculation 
of the kSZ effect we assume that the nonlinear density 
ffuctuations are un corrected with the b ulk velocity field 
in which they lie. Zhang et al. (2004) argue that this 
approximation may not hold m highly nonlinear regimes 
where contributions from the curl of the nonlinear veloc- 
ity field may become important however we neglect these 
corrections here. 

As in the previous section, we find a power-law approx- 
imation for estimating the kSZ power as a function of cos- 
mological parameters. The kSZ angular power spectrum 
was calculated under the full analytic formulation for a 
large suite of WMAP7-allowed LCDM cosmologies. An 
MCMC was then performed in the 6 dimensional fitting- 
function parameter space and best-fit marginalized val- 
ues were found and are listed in Table [T] 

In Figure |3] we compare our calculation of the kSZ 
power spectrum (solid red line) with that measured from 
recent simulations (black lines). As in Figure [2] we 
plot the power spectrum predicted by the Ma re-Nostrum 



simula tion (dotted), the 'standard' model of Trac et al. 
(2010) (dot-dashed) and the non-radiative hydrodynam- 



ical simulations of Battaglia et al. (2010!) (solid). The 
red dashed line shows our model tor the contribution to 
the kSZ signal from inhomogeneous reionization (not in- 
cluded in the other lines), as described in the following 
section. 

2.3.2. Patchy Reionization 

We also consider the contribut ion to the kSZ power 
from inhom ogeneous reionization (Gruzinov & Hu "1998^ 
, Knox et al. 1998; Hu 2000; Zahn et al. 2005; McQuinn, 
et al. 2005; Iliev et al. 2007j). Simulations and analytic 
models of HII bubble formation both indicate that the 
first galaxies and quasars were highly clustered and led 
to gradual reionizatio n in "bubbles" th at quickly grew to 
sizes of several Mpc ( Zahn et al.||2010 ). 



TABLE 1 

SZ COSMOLOGICAL SCALING 



i 


Anv fuK^l 

^^<JV Lr^ J 


Tig 


f2j, 




^ o 


7" 






0.96 


0.045 


0.22 


0.8 


0.09 


500 


1.18 


-1.44 


1.83 


-1.06 


4.36 


0.25 


1000 


1.81 


-1.36 


1.91 


-1.13 


4.82 


0.24 


2000 


2.64 


-0.94 


1.96 


-1.12 


5.26 


0.22 


3000 


3.06 


-0.45 


1.94 


-1.03 


5.38 


0.20 


4000 


3.33 


-0.22 


1.96 


-1.04 


5.54 


0.17 


5000 


3.53 


-0.03 


1.98 


-1.05 


5.66 


0.15 


6000 


3.67 


0.13 


2.00 


-1.06 


5.75 


0.13 


7000 


3.78 


0.27 


2.01 


-1.07 


5.83 


0.12 


8000 


3.87 


0.38 


2.02 


-1.08 


5.89 


0.10 


9000 


3.94 


0.49 


2.03 


-1.09 


5.94 


0.09 


10000 


4.00 


0.58 


2.04 


-1.10 


5.99 


0.08 


£ 


^tbZ; L/-^^^ J 






n 




h 






0.96 


0.045 


0.265 


0.8 


0.71 


500 


2.06 


-1.01 


2.41 


0.69 


8.57 


1.30 


1000 


3.59 


-0.75 


2.45 


0.63 


8.49 


1.43 


2000 


4.86 


-0.36 


2.52 


0.53 


8.40 


1.65 


3000 


5.04 


-0.08 


2.57 


0.48 


8.36 


1.73 


4000 


4.82 


0.15 


2.62 


0.44 


8.34 


1.88 


5000 


4.44 


0.31 


2.66 


0.42 


8.33 


2.03 


6000 


4.03 


0.49 


2.70 


0.39 


8.32 


2.13 


7000 


3.62 


0.58 


2.73 


0.38 


8.32 


2.18 


8000 


3.24 


0.77 


2.78 


0.36 


8.32 


2.27 


9000 


2.89 


0.87 


2.80 


0.35 


8.33 


2.32 


10000 


2.59 


0.96 


2.83 


0.34 


8.33 


2.37 



Note. — The ^-dependent power-law cosmological scalings for 
the Ostriker- Vishniac effect and the thermal SZ effect. The num- 
bers immediately below the cosmological parameters are the pivot 
points for the power law, and the numbers in the table are the 
power-law indices. For example, the top row says that for the OV 
effect, D500 = 1.18mK2 {ns/0.96)~^-'^'^{nt/0m5y-^'^... 



In 



our estimat es we use the an alytic Monte-Carlo 
"FFRT" model oflZahn et all (120101). It has been shown 



to agree well with the most sophisticated radiative trans- 
fer simulations on scales of 100 comoving Mpc/h, while 
having the added advantage of allowing the modeling of 
arbitrarily large volumes (the analytic scheme is about 
4 orders of magnitude faster, at a given dynamic range, 
than radiative transfer). This is especially important for 
kSZ, since large scale velocity streams lead to the bulk 
of the signal. Our particular template (dashed red line 
in Figure 3]) was calculated in a 1.5 Gpc/h cosmological 
volume where x- and y-axes correspond to roughly 15 de- 
grees on a side and z-axis corresponds to redshift, with 
a median redshift of 8. We shift this template left-right 
logarithmically by a "patchy shift" parameter Rp; that 
is, the power spectrum for a given shift is related to the 
fiducial Rp — 1 spectrum by, 

Cf, = C 



fid 

Rpxe 



(13) 



Rp is to be thought of as scaling the size of the bubbles, 
and is, to good approximation, proportional to the du- 
ration of the patchy phase. The timing of reionization 
has a secondary small effect on the shape and amplitude, 
which we neglect here. 

2.4. tSZ-DSFG Correlation 

It is reasonable to expect some correlation between the 
DSFG clustering and tSZ components since they both 
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trace the same underlying dark matter distribution, with 
significant overlap in redshift. Simulations which asso- 
ciate emission with individual cluster member galaxies 
predict anti-correlations — DSFGs fill in SZ decrements at 
frequencies below 217 GHz considered here — on the order 
of tens of percent, w ith a correlation coe fficient nearly in- 
dependent of scale ( Sehgal et al. 2010). This effect was 
explored in yhirokoif et al. 12010) which found correla- 
tion consistent with zero but with significant uncertainty 
due to degeneracies with the tSZ and kSZ components. 

Assuming a fixed correlation rtgz.c^ the total power 
spectrum isn't simply the sum of the tSZ and DSFG 
clustering terms, but must also include a term given by, 

C,,... = nsz,c [^JclJWu' + y^S^^) (14) 

Note that this effect can be larger than either compo- 
nent individually in cross spectra between frequencies in 
which each component is large. For example, even with 
only moderate levels of correlation, the 217 x 70 GHz 
foreground contribution would be dominated by this cor- 
rection because of the large DSFG power at 217 GHz and 
tSZ power at 70 GHz. 

2.5. CMB 

For the primary CMB signal itself, we use an 8- 
parameter model which includes the baryon density 
r2f,/i^, the density of cold dark matter ilch^, the optical 
depth to recombination r, the angular size of the sound 
horizon at last scattering Q, the amplitude of the pri- 
mordial density fluctuations ln[10^°As], the scalar spec- 
tral index Us, the dark energy equation of state parame- 
ter w, and the tensor-to-scalar ratio r. Freeing w opens 
up the "geometric degeneracy" which is typically broken 
by adding an external dataset such as supernovae data. 
Rather than do this, we simply put a ±0.3 Gaussian prior 
on w to reasonably constrain the chain while allowing it 
to explore the parameter space. 

Because we are interested in the simplest description 
of how the foregrounds affect cosmological parameters, 
we do not consider extensions to our model such as a 
running spectral index, non-fiat universes, non-standard 
effective number of neutrino species, or a difference in 
primordial helium from standard big bang nucleosynthe- 
sis. Due to the small angular scales where they affect 
the CMB anisotropy, it is possible that such parameters 
could be even more degenerate with the foregrounds than 
the "vanilla" set we consider. 

For quick and highly accurate CMB calculations dur- 
ing ou r MCMC chains, we use a PICO ( |Fendt fc Wandelt 
2007| interpo lation of a training set generated by CAMB 
( [Lewis et al.]|2000l ). Though more recent recombination 
codes exist, we expect the one used in our training set to 
be sufficient for our forecasting. Additionally, the train- 
ing set includes the option of a non-linear lens ing contri- 
bution described in Challinor & Lewis (2005) which we 



3. FIDUCIAL MODEL AND CURRENT CONSTRAINTS 

For our forecasting, we create simulated power spec- 
tra (henceforth the "simulated data") using the model 
described in the previous sections. We pick one single 
set of model parameter values, called the "fiducial val- 
ues" or the "fiducial model" in general, which will is the 



baseline for the different cases of simulated data which 
we consider. The model used to analyze the simulated 
data, which generally contains small changes relative to 
the fiducial model, will be called the "analysis model." 

In Table [3] we summarize all of the parameters in our 
model, the naming convention, and their fiducial val- 
ues. The fiducial values are chosen to be consistent with 



current cos mological constrains from WMAP7 (Komatsu 
et al.|2010 ), and with constraints on the foreground com 



ponents from ground-based data such as SPT and ACT. 
In the following paragraphs, we describe the method used 
to arrive at our fiducial model. 

Because the expected SZ power depends strongly on 
cosmology, special care was taken so that our fiducial 
SZ power and cosm ology agree. To achie ve this, we use 



the constraint from Lueker et al. (2010) 
combination Vtsz 



on the linear 

4:t±i.5 AiK2 (at e : 



.46 X Vksz 

3000 and 153 GHz) along with the cosmological scalings 
in Table [T] We then importance sample the WMAP7 
ACDM-^TENd^ chain by calculating at each step the 
expected SZ (kmetic and thermal) power assuming no 
theory unce rtainty, then applying the prior from Luekerj 
et al. ( 2010[) . The new best fit point in the post- processed 
chain mainly shifts SZ power up relative to best fit SPT 
value and as down relative to the best fit WMAP7 value. 
All other cosmological parameters are also affected (at 
a smaller level), and their new mean values form our 
fiducial cosmology, which remains 1-cr consistent across 
all parameters with WMAP7. 

For the radio sources, the tightest constraints on the 
expected Planck power come from the Vieira et al. ( 2010 1 
catalog which contains sources in the decade o t bright- 
ness just be low the Planck flux cut. Fitting a |de Zotfr 



et al. ([2005| model to the data yields the values listed 
in Table [3j notably radio Poisson power of 133 /iK^ at 
143 GHz assuming a 330 mJy flux cut. 

Since the DSFG Poisson contribution is nearly inde- 
pendent of flux cut, we expect the same Poisson power in 
Planck maps as in SPT maps, adjusting only for band- 
pass differences. We get our fiducial value for Planck 
Poisson power at 143 GHz by extrapolating in frequency 
from th e best-fit value of the SPT 150 GHz power as 
given in iShirokoff et al.| (|2010|) . Our fiducial values for 



ap and ag also come ifrom the best-fit values in |Shi-| 
' ' " ' ~ W e set a p to 0.4 following the argu- 

(20041, although it is not yet well 



rokoff et al. (^2010 



ments in Knox et al 
constrained by observations. We adopt a clustering tilt 
uc — I so that it (roughly) has the shape expected at 
small scales due to the the observed clust ering properties 
of hig h-redshift galaxies. As argued by Scott fc White 
(11999!), the observed clustering properties of z ~ 3 Ly- 
man break galaxies, namel y an angular correlatio n func- 
tion proportional to 9~^'^ ( Giavalisco et al.|[i998 ), corre- 
spond to Di (x £^'^. Since we multiply the power-law by 
the linear theory template, this is similar (at i > 1500) 
to the power-l a.w only oc sha pe u sed as baseline 



models in both Dunkley et al. (2010) and Shirokoff et al.| 
( 2010 1 . Although ruled out by the Planck data, our fidu- 
cial niodel is at least closer to the measurements than all 
of the other models plotted in Figure[Tl The agreement is 
sufficient for our purposes here, thoiigh we will certainly 



Available at http://lambda.gsfc.nasa.gov/ 
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TABLE 2 
Survey Properties 



Band 


T CF-/R~I 
1 (rj/D) 


Bccim 




(GHz) 


III l\ yr'vvi ivi\ 


1 Oil l^llilll 1 




Planck 








1 u 




14 


f , — 70% 
J sky — 'U/o 


100 


61 (98) 


10 


e _ oon ^ T 

Oc — oo\j Liii) y 


143 


42 (80) 


7.1 




217 


64 (132) 


5 




Ground- Deep 








90 


53 


1.6 


f,ky = 100dcg2 


150 


13 


1.15 


Sc = 6.4 mjy 


220 


35 


1.05 




Ground- Wide 








90 


53 


1.6 


fsky = 1000 deg2 


150 


18 


1.15 


Sc = 6.4 mJy 


220 


80 


1.05 





Note. — Instrument properties used to generate simulated 
power spectra. The beam width is given as a fuU-width-half-max. 
Sc refers to the flux cut above which brighter sources are masked 
out. 

be updating our CIB mode l ing in the near future 



Following Battye et al. 



(2010) which found a mean 
fractional polarization ot 4.5% at 86 GHz (and varying 
weakl y with frequency) fo r the WMAP point source cat- 
alog (Wright et al. 20091, we adopt a fiducial value of 
fn = 0.05. b'or USFGs we expect an even smaller level 
of average polarization fraction. Polarized dust emission 
arises due to alignment of grains in interstellar magnetic 
fields. We somewhat arbitrarily set fo — 0.01 for our 
fiducial model which is consistent with the finding that, 
in our own galaxy, the coherence length for magnetic 
fields is much smaller than the extent of the dust emis- 



sion fPr unet et allll998| 

Figure |4] shows the haucial CMB and foreground con- 
tribution to Planck TT, TE, and EE power spectra (with 
the exception of tSZ-DSFG correlation which is plotted 
at 30% rather than its fiducial value of 0%). 

4. SURVEY PROPERTIES 

We consider simulated Planck data in the four bands 
between 70 GHz and 217 GHz. These are chosen because 
they contain nearly all of the significant CMB informa- 
tion. Though the neglected channels place little extra 
constraints on the CMB, they are crucial for understand- 
ing and cleaning the foregrounds. We consider their ef- 
fect implicitly by testing limits such as lowered Poisson 
power amplitudes, or fixed DSFG clustering shapes. Ad- 
ditionally, we also consider the benefit of higher reso- 
lution ground-based data, which we model after SPT 
90 GHz, 150 GHz, and 220 GHz channels. We divide 
the data into two fields: a 100 deg^ "deep" field and a 
1000 deg^ "wide" field. Wc henceforth refer to these two 
datasets as Ground-deep and Ground-wide. The depths, 
sky coverage, and flux cuts used in our forecasting are 
summarized in Table [2] 

Our simulated data take the form of auto and cross 
spectra from as many bands as are present for a given 
patch of sky. The four Planck frequency channels form 
10 TT, EE, and BB and 16 TE power spectra, with an 
additional 18 TT power spectra from the three extra fre- 
quencies in regions of Ground overlapp] We do not as- 

In general N frequency channels can be used to create N{N + 



sume overlap between Ground deep and wide, nor do 
we form cross spectra between Ground temperature and 
Planck polarization as these are expected to be a very 
small contribution to the CMB and foreground informa- 
tion. Planck BB polarization is also ignored except in 
one test case where we find its impact is minimal on our 
cosmological parameterization. 

We simulate power spectrum assuming a uniform 
masking threshold across the sky. The only exception 
is in the case of Planck and Ground overlap. For such 
patches of sky, we assume Planck maps can be masked 
using a point source mask from the higher resolution 
Ground data. Thus, for the overlap areas, even the 
Planck auto spectra will have greatly reduced radio Pois- 
son power. 

The non-zero width of frequency bandpasses creates a 
different effective frequency for each component in each 
band. For components with uncertain spectral shapes, 
the variation in effective frequency leads to percent level 
corrections which can be neglected. In this paper, val- 
ues quoted from ground-based experiments are normal- 
ized at, and explicitly cite, the corresponding effective 
frequency. For Planck, it is sufficient for our forecasting 
purposes to ignore this and use the nominal band centers 
for all components. 

5. FORECASTING METHODOLOGY 

The analysis of the simulated data assumes per- 
fectly known Gaussian beams, no calibration uncertainty, 
isotropic noise, and ignores the effects of mode- mode cou- 
pling on the cut sky. While these assumptions are not 
sufficient for modeling real data, we expect them to be 
adequate for our purpose of modeling the extragalactic 
foregrounds, and understanding their importance on cos- 
mological parameter biases and statistical errors. 

Under these assumptions, the so called "pseudo power 
spectrum" which includes both signal and noise is given 

by, 



(15) 



Here i and j each label one of the maps, and the noise is 
parameterized by the weight per solid angle w, and the 
beam width in radians ai,. To calculate the error on our 
estimates of the signal power spectrum, we use a Fisher 
matrix approximation which is nearly exact in the £ range 
we consider. The resulting "bandpower covariance" is 



(2/ + l)fsky 



{CilCjk + CikCji) (16) 



where we have suppressed the £ dependence for nota- 
tional simplicity. We use the covariance to form the like- 
lihood as a function of parameters 9, 



2 In £(61) 



(0) 



'^kl 



Note we have neglected the normalization term since it 
does not vary with 9. 

Our simulated data are the mean expected power spec- 
tra; i.e., they do not include a sample of the errors from 

l)/2 power spectra of type TT, EE, and BB, and of type TE 
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Fig. 4. — All 36 power spectra which can be formed from Planck 70GHz-217GHz temperature and E-mode polarization, and the 
prediction of our fiducial model for the CMB and foreground power in each of them (with the exception of the tSZ-DSFG correlation which 
is shown at 30% instead of its fiducial value of 0%). The black dashed line shows the errors bars for £-h'ms of width A£ = 256. Dotted lines 
indicate negative power. 
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TABLE 3 
Summary of Model Parameters 



Parameter 


Fiducial Value 


Current Constraints (1(t) 


Definition 


Cosmological 


.022565 


.00073 


Baryon density 




.10709 


.0063 


Cold dark matter density 


e 


.010376 


.000029 


Angular size of the sound horizon at last scattering 


r 


.0799 


.015 


Optical depth to reionization 


w 


-1 


.13 


Dark energy equation of state parameter 


rig 


.9669 


.014 


Scalar spectral index 


ln(10^''As) 


3.1462 


.045 


Scalar amplitude 


T 


. io 




Tensor-to-scalar ratio 


Dusty Poisson 








ao 


3.8 


0.35 


Spectral index 


<^D 


.4 




Spectral index intrinsic spread 


■Dd 


5.9 /iK^ 


0.8 


Amplitude at ^ = 3000, u = 143 GHz 


in 


.01 




Dusty polarization fraction 


Radio Poisson 










-.5 


0.1 


Spectral index 




.1 


< 0.6 (95%) 


Spectral index intrinsic spread 




133 


27/iK2 


Amplitude at f = 3000, v = 143 GHz, 5c = 330 mjy 


IR 


-.8 


0.1 


Brightness function power law index 


Ir 


.05 




Polarization fraction 


Dusty Clustered 








ac 


3.8 


0.4 


Spectral index 


-De 


3.9 AtK^ 


1.2AtK2 


Amplitude at ^ = 3000, u = 143 GHz 


nc 


1 




Nonlinear tilt 


SZ Effects 








'Dtsz 


4.3 AtK^ 


< 6.8 /iK2 (95%) 


tSZ amplitude at ^ = 3000, v = 143 GHz 


'DkSZ,OV 


2.7aiK2 


< 6.5AiK2 (95%) 


OV amplitude at ^ = 3000 


'DhSZ,p 


1.5 AtK^ 


< 6.5AiK2 (95%) 


Patchy amplitude at £ = 3000 


Rp 


1 




Patchy shift 


Correlations 








rtSZ,C 







Correlation between tSZ and DSFGs at £ = 3000 



Note. — A summary of the parameters in our model. The fiducial values generate our simulated data. The current 
constraints column gives the la constraints on our model given WMAP power spectra and radio source counts, SPT 
power spectra and radio/DSFG source counts, and ACT power spectra. Note that due to the process by which the 
fiducial values were chosen (Sec.jsl they are not necessarily the most likely values given current data; they are, however, 
totally consistent with the most likely value to within la. 
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the bandpower covariance matrix. Leaving out these 
fluctuations has the benefit of making the best fit 
equal to exactly (as long as our analysis model and 
simulation model are the same) and has no affect on our 
forecasting abilities. 

Tests we performed showed that a Gaussian propaga- 
tion of uncertainty from the Cis to the model parame- 
ters can be insufhciently accurate due to the highly non- 
Gaussian shape of the foreground parameter posterior 
likelihoods. Instead, we run a full Markov-Chain Monte- 
Carlo (MCMC) analysis using a custom mult i-frequency 
extension to CosmoMC ( jLew'is &: lBridlel2002| ). The code 
for this extension along with chains tor the results quoted 
in this paper are available onlinfl^ 



6. COMPRESSION TO A CMB POWER SPECTRUM 
ESTIMATE 

Before getting to our results, it is useful to explore the 
foreground contamination in a more model-independent 
manner, motivated by two drawbacks of our procedure. 
First, there is a large amount of data one must work 
with — our bandpower covariance matrix at each £ is 
46x46 for Planck and 64x64 for Planck+Ground. Sec- 
ond, 17 foreground parameters must be marginalized 
over, and if one wanted to examine constraints on a new 
cosmological model, the whole procedure would have to 
be repeated. Here we present a procedure for compress- 
ing all the power spectra to 1) a single CMB estimate 
and 2) a low dimensional parametrization of the residual 
foregrounds in this estimate. We describe the procedure 
here for temperature-only power spectra with errors that 
are uncorrelated from multipole to multipole. The gen- 
eralization to include £ to £ correlations and polarization 
is in Appendix [Xj 

Given N power spectra (for example the 10 TT spec- 
tra we consider for Planck) ^ we would like to split our 
data up into A'^ — 1 linear combinations of power spec- 
tra that have no sensitivity to the CMB and then find 
the remaining linear combination that contains CMB and 
whose errors are uncorrelated with those of the N — 1. 
With this split made, we can then derive our foreground 
constraints using the CMB-free linear combinations. Do- 
ing so means foreground constraints can be made inde- 
pendent of our modeling of the CMB (other than the 
assumed frequency dependence). 

We use the N — 1 CMB-free linear combinations to 
find the constraints they place on our foreground model 
parameters via MCMC. For each point in the chain we 
can determine the contribution to the CMB linear com- 
bination. We sample over all these contributions to find 
the mean contribution and fluctuations about that mean. 
We find a low-dimensional description of the fluctuations 
via a principal component decomposition. 

6.1. Splitting the power spectra into CMB-free and a 
CMB estimate 

Let us begin by first considering arbitrary linear com- 
bination of the power spectra, 

C"^ = ^<C, (18) 

i 

http : // student . physics . ucdavis.edu /^millea/data /millea20 1 1 




1000 1500 2000 2500 3000 3500 

; 

Fig. 5. — Top: The ^-d epe ndent weightings which form the CMB 
linear combination (Eq. |28| l. All possible auto/cross spectra from 
Planck channels in Table were considered. Dashed lines indicate 
negative weight. Middle: The mean foreground contribution to the 
CMB linear combination for our fiducial model. Note for example 
that tSZ (purple) is not present at high £ because only 217 GHz is 
used there. The dominant non-Poisson component for the l-range 
where Planck is most sensitive is the DFSG clustering. Bottom: 
Principal components of foreground residuals (constrained by the 
CMB-free linear combinations) with amplitudes set to 1-a. Note 
that we only need two principal component amplitudes to be ac- 
curate to > l^iK^. (The errors in bin widths of Ai = 256 for both 
the CMB linear combination and for 217 GHz alone are plotted as 
dashed lines in the bottom two plots.) 

where Ci are the i = 1 to N power spectra. The weight- 
ings we consider will be ^-dependent; the lack of any la- 
beling by £ is solely for notational simplicity. Here, /i is 
merely a label to distinguish different weightings; the 
are a linear combination of the old power spectra with 
weig ht wf . Note that if a weighting satisfies. 







(19) 



it is not sensitive to the CMB. 

We would first hke to find the CMB weighting ui^mb 
which will be statistically orthogonal to the — 1 linear 
combinations which satisfy the CMB-free condition (Eq. 
19). We would also like this weighting to be properly 
normalized so that. 



1. 



(20) 
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To satisfy the orthogonahty condition it helps to work 
in a primed space defined by a hnear transformation via 



w 



!W7 



(21) 



where L is the Cholesky decomposition of the bandpower 
error covariance matrix, E = LL^ . The advantage of 
the primed space is that the basis vectors in the primed 
space correspond to power spectra whose errors are sta- 
tistically orthogonal; i.e., with the weightings set so that 
S/^ia (now setting fj, = 1...N) the errors in the 



w, 



corresponding power spectra satisfy 

{SC'^SC") = 5^,. (22) 

The primed weights that satisfy the CMB-free condi- 
tion satisfy 



Thus any power spectrum with primed weighting pro- 
portional to 



is perpendicular to any vector satisfying the CMB-free 
condition, as one can easily verify. To find the CMB 
weighting in the unprimed space we perform the inverse 
transform and no rmalize to satisfy the normalization 
condition (Eq. [20]) 



,,,CMB _ 



(25) 
(26) 



Note that this is the expression for inverse-variance 
weighting. 

Our remaining task is to construct the TV — 1 CMB-free 
weightings in a manner that leaves them all statistically 
orthogonal to the CMB weighting. We do so by apply- 
ing the Gramm-Schmidt procedure in the primed space. 
This gives us iV — 1 orthogonal vectors that are all or- 
thogonal to the CMB direction as well, that we will call 
v'J^ for ^ — 2,N. The weightings in the unprimed space 
are then given by 



v ^ 

I la a 

We now define the matrix W so that 



: W. 



CMB 



W^,^vt{foi ^i = 2...N}. 



(27) 



(28) 



This matrix defines the linear combinations of the power 
spectra that have all the properties we desire. The first 
row is the optimal CMB weighting and subsequent rows 
give the — 1 CMB-free linear combinations. All the 
linear combinations are statistically orthogonal; i.e., the 
covariance matrix for the new power spectra, W'^T,W, 
is diagonal. Furthermore, W is non-singular so we have 
not lost any information through this re-weighting. 



6.2. Modeling the foreground residuals with principal 
components 

With the weight matrix W defined, we can constrain 
the foreground power in the A^ — 1 CMB-free power spec- 
tra by running an MCMC chain. Despite the large num- 
ber of parameters and power spectra, this analysis is fast 
in practice because the foreground model consists of sim- 
ple analytic forms and precomputed templates, and does 
not depend on any costly Einstein-Boltzmann solver or 
lensing models. For the set of foreground parameters 
at each step in this chain, we calculate the correspond- 
ing foreground contribution to the CMB linear combina- 
tion. These ^-dependent contributions form the columns 
of the Y matrix in a principal component analysis (see 
Appendix [B I. Following the PC A procedure, we have a 
few principal components and priors on their amplitudes 
which must be marginalized over in a separate chain 
which uses only the CMB linear combination. 

6.3. Discussion of linear combination analysis 

The weightings which make up the CMB and CMB- 
free linear combinations depend on the bandpower co- 
variance matrix, and thus on the noise properties of the 
instrument, any filtering which is performed, and on the 
true power spectrum on the sky. The principal compo- 
nents for the foreground residuals also depend on the 
choice of foreground model. For a Planck temperature- 
only forecast and for our fiducial model, we present the 
results of a linear combination analysis. 

In the top panel of Figure [s] we plot the weights for the 
CMB linear combination as a function of £. At high £ 
where the measurement is noise dominated, nearly all 
of the CMB information is contained in the 217 GHz 
map which is the least noisy. At lower £ where we be- 
come dominated by cosmic variance, the CMB informa- 
tion comes from the channels with the lowest foreground 
contamination. 

Given these weights, we plot in the middle panel of 
Figure [5] the foreground contribution to the CMB linear 
combination and the error bars on this new power spec- 
trum. Also shown are the error bars for the 217 GHz 
channel alone for comparison. The maximum improve- 
ment is at £ = 2000 where the error bars tighten by a 
factor of 1.4. We also see that the dominant contribution 
to the foreground power in the £ range where Planck is 
most sensitive to the CMB is the radio Poisson, followed 
by the DSFG clustering. 

Finally, we perform a PGA on the foreground residuals 
in the CMB linear combination. The first several princi- 
pal components are shown in the bottom panel of Figure 
[sj We find that all of the variation > 1 /j,K^ can be de- 
scribed by two amplitude parameters, as compared to the 
14 parameters which govern these foregrounds. Another 
way to put this is that using the CMB-free combina- 
tions we can clean out almost 140 /iK^ of foregrounds (at 
£ = 3000), leaving only tens of /iK^ of residual uncer- 
tainty, modeled with the two principal components. 

7. RESULTS 

With the model and forecasting tools in place, we are 
ready to present the results of our main analysis. We 
want to find which components can potentially cause 
large biases in an analysis of Planck data, so that we 
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can model them with sufRcient care. We would also like 
to know how much constraining power is reduced due to 
foreground confusion. Could significant improvements in 
cosmological parameter constraints be achieved by using 
additional data and/or modeling? To answer these ques- 
tions, we run a suite of forecasting analyses aimed at 
singling out the effects of each foreground contribution. 

The next sub-sections are organized as follows. First 
we examine the importance of each component by turn- 
ing it on or off in the analysis. For the components which 
we deem important, we check whether our modeling is 
sufficient to protect the cosmological parameters from bi- 
ases, both at the Planck and PZancA; -|-Ground sensitivity 
levels. We then examine the degradation in statistical er- 
rors from the need to marginalize over the foregrounds, 
and finally we explore the impact of adding in ground- 
based data. 

7.1. Importance of the Different Components 

Figure |6] shows the effect of removing four foreground 
components — the DSFG clustering, tSZ, kSZ, and tSZ- 
DSFG correlation — one at a time from the analysis 
model, while they are actually present in the simulated 
data at their fiducial value. We present the results by 
plotting likelihood contours in the ris and flch^ plane, 
since changes in those two parameters affect the primary 
CMB at the smallest scales and are the most suscepti- 
ble to foreground biases. We also show the amplitudes 
of the clustering and SZ effects as their ^-space shapes 
make them most degenerate with cosmological parame- 
ters. All of the chains in this section include only Planck 
power spectra in the simulated data. 

We expect the DFSG clustering to be extremely im- 
portant to model since it is the second largest fore- 
ground contribution to the CMB linear combination 
in the €-range where Planck is most sensitive. When 
marginalized over, this contribution is constrained to be 
10.5 ± 0.6 yuK^ at i = 1500, so setting it to zero is about 
an 18(7 change. The dot-dashed green contours in Figure 
[6] show that this is compensated by a systematic bias of 
7(7 is Us and llcr in flch^, along with an increased kSZ 
power to about 30 fiK^ . Using the middle panel of Figure 
[5] as a visual guide, we can examine how this happens. 
Though the kSZ increases to compensate for the miss- 
ing clustering power at high-^, its shape is flatter than 
the DSFG contribution to the CMB linear combination, 
so Us decreases to roughly remove the extra power at a 
low-£. 

At 217 Ghz, the kSZ power in our fiducial model is 
a factor of 20 times smaller than the DSFG clustering 
and is therefore (in terms of cosmological parameter es- 
timation) less troublesome. However, due to its identical 
frequency dependence to the CMB, we do expect the am- 
plitude of the kSZ signal to be degenerate with cosmo- 
logical parameters. The solid black contours in Figure 
[6] demonstrate the bias introduced if the kSZ compo- 
nent is omitted from the foreground modeling. We find 
a roughly 0.5cr bias in Ug as it increases to try to fill in 
the missing 4.2 fiK^ of kSZ power. The DSFG and tSZ 
plane (bottom right panel of Figure [6| shows that the 
other foregrounds are largely unaffected. 

The thermal SZ component is neither frequency inde- 
pendent, nor does it contribute as much power as the 
DSFG clustering, so we do not expect a bias as large as 
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Fig. 6. — 68% (and 95% in the bottom panel) confidence con- 
tours for a suite of test cases examining the effect of neglecting to 
model different foregrounds. Unless explicitly stated above, other 
parameters were included in the data at their fiducial values listed 
m Table [3] and were marginalized over in the analysis. A''^ corre- 
sponds to the maximum number of £-bins per power spectrum one 
coul d us e and still detect the error in modeling at 3 sigma (see 
Sec. |7.l| for further discussion). 

in either previous case. It does, however, project into the 
CMB linear combination to an £-shape very similar to the 
CMB itself (see the middle panel of Figure [5]), making it 
more likely to be degenerate with cosmological parame- 
ters. From the results, we see about O.Sfi biases in each 
of Us, ^ch^, and ilbh^. 

Finally, we consider neglecting a 30% tSZ-DSFG cor- 
relation, a value on the high end of expected correlation. 
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Fig. 7. — The effect on cosmological parameters from tr ying to fit 
our m odel to simulated data which includes ( orange) the|Battaglia 
let al. f2010 ) tSZ template and (black) the [Amblard L Uoor^ 
l |l^UUY I clustering template. These two models are tue most dissim- 
ilar to ours, and thus show our model can protect against biases 
of a few percent up to Planck sensitivity. The inclusion of Ground 
data necessitates more detailed modeling of only the clustering. 



but still consistent with [Shirokoff et al.| ( |2010[ ). We ex- 
pect this to have the smallest ettect on the cosmological 
parameters since the power contribution is sub-dominant 
to all of the other foreground components at all fre- 
quencies which appear in the CMB linear combination 
at > 1%. While the measured tSZ amplitude is biased 
at a few sigma as it raises to compensate for the miss- 
ing power, the effect is not large enough to significantly 
impact any of the cosmological parameters. 

One question is whether any of these analysis errors 
would be caught by a goodness-of-fit test. To address 
this question we present Ax^ values in the table in Fig. [jj 
We can expect rms fluctuations in to be y^Ni, where 
Ni, is the total number of bandpowers which is roughly 
equal to the number of degrees of freedom. If one is 
searching for signs of a contaminant that is very slowly 
varying in £, then one would bin coarsely to reduce the 
statistical fluctuations in x^, to make a more stringent 
goodness-of-fit test. 

We define Ng to be the number of £-space bins such 
that the absolute Ax^ from the fit is 99.7% inconsis- 
tent with random fluctuations. Thus we have N{ — 
(Ax^)^/(9 X A^spoc) where A'spcc is the number of power 
spectra (36 here). We see that binning would not have 
to be coarse at all to detect the poor fit caused by ne- 
glecting clustering. We also see that for the other entries 
in the table, binning would have to be extremely coarse 
for the fits to be noticeably poor. Indeed, the binning 
would have to be coarser than is practical since the sig- 
nals of interest, as well as the contaminants, would vary 
significantly across a bin. We conclude that only the "no 
clustering" case would produce a noticeably bad fit for 
Planck only. 

7.2. Modeling Sufficiency 

Given the demonstrated importance of the foreground 
components, we would now like to see if our modeling 
is sufficient to protect the cosmological parameters from 



biases if we have modeled the components, but modeled 
them incorrectly. In this section we consider the DSFG 
clustering and the tSZ effect. 

For the DSFG clustering, we turn to the models plot- 
ted in Figure [T] Our parameterization s hould have the 
most trouble reproducing the Amblard & Cooray ( 2007 1 



model, which switches to a power-law (as a consequence 
of non-linear clustering) at £ ~ 2500 rather than at 
£ = 1500 as in our fiducial model. The orange contours 
in Figure [7] show the results obtained when fittin g our 
model to simulated power spectra that assume the Am-| 
[blard fc Cooray| ( |2007[ ) clustering template. As we had 
hoped, tor the case ot Planck only (solid lines), there is 
no significant biasing. 
We also explore the ability of our tSZ princ ipal compo 



nent model (based on the analytic model of Shaw et al. 
20101 to encompass the variations in the tSZ models 



shown in Figure^ We elect the Battaglia et al. (20101 
one as the most dissimilar, since it lacks the ettects of ra- 
diative cooling, an d should be the most difficult for the 
Shaw et al. ( 2010[ ) model to reproduce. Despite these 
ditterences, l^'igure [T] shows that for Planck the model is 
sufficient to encompass the shape uncertainty and pro- 
tect cosmological parameters. 

When we add in Ground (dashed lines), the require- 
ments on the modeling accuracy are more stringent. For 
the clustering case, we see an almost Icr bias from using 
our fiducial mo del when the true model is the Amblard] 
& Cooray (j2007) clustering template. Analyses with cur- 



rent da ta can tolerate muc h more discrepant clustering 
shapes (Dunkley et al. 2010 ). For future P/ancfc -I- Ground 
analyses, the clustering shape will need to be modeled 
more accurately. For tSZ the modeling appears to be 
more robust; tSZ-induced biases are small even in the 
Planck+Gromid case. 

7.3. Statistical Error Increase with and without 
Auxiliary Data 

We have demonstrated the possibility of cr-level bi- 
ases in cosmological parameters arising from failure to 
model foregrounds. To prevent these biases, the fore- 
grounds need to be jointly estimated or marginalized 
over. We now turn to two questions: 1) How much do 
the cosmological parameter statistical errors degrade due 
to foreground uncertainty? and 2) How much can be 
gained from using other data to constrain foregrounds 
and thereby reduce that degradation? 

The top two rows of Table |4] show the effect of 
marginalizing over our entire foreground model as op- 
posed to fixing it at fiducial values. In each row, the 
difference from 100% is the percent degradation due to 
foreground marginalization. The second row shows the 
degradation is limited to 20% for Ug and 10% for As and 
rt/jh^. We see no degradation in r and r since they are 
mainly constrained by large scales where the extragalac- 
tic foregrounds we consider are negligible. The dark en- 
ergy equation of state w is unaffected because it is mainly 
constrained by our ±0.3 prior. 

Ground data can help reduce this degradation by bet- 
ter constraining the foregrounds using auto and cross 
spectra that are more sensitive at small scales. The im- 
provement from adding these to the simulated data is 
shown in the row labeled PZancfc-f Ground. The mea- 
surement of DSFG shot noise is improved ten-fold, with 
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TABLE 4 
Statistical Error Dec;radation 












ris 








»C 




T^kSZ 


Planck (fgs fixed) 


1.1x10--* 


1.0x10-=* 


2.6x10-4 


3.0x10-3 


1.3x10-2 












Planck (fgs marginaUzed) 


110% 


100% 


100% 


120% 


110% 


3.4 


6.0 


1.3 


1.0 


4.4 


P/ancfc+Ground 


100% 


100% 


100% 


110% 


100% 


10% 


50% 


40% 


60% 


60% 


Planck (Clean DSFG) 


100% 


100% 


100% 


110% 


100% 


60% 


100% 




100% 


100% 



Note. — Entries in bold show the standard deviation for each parameter for that particular row. Subsequent percentages reference 
the bold entries in the same column. Dashes indicate the parameter was fixed, while blanks mean the parameter is not applicable to that 
case. The normalization parameters T>x are in units of /iK^. The different cases correspond to: (fgs fixed) Fixing all of the foregrounds at 
their fiducial values, (fgs marginalized) Marginalization over our full foreground model. (+Ground) Also including Ground auto and cross 
spectra in the simulated data. (Clean DSFG) Assuming 90% reduced clustering power due to cleaning from higher frequencies. 



the clustering and SZ effects also tightened by a factor 
of two. The radio amplitude is improved through con- 
straints on the spectral dependence, and could be fur- 
ther improved though a prior on 7^ from Ground source 
counts. The effect on the cosmological parameters is to 
remove essentially all of the degradation we incurred from 
marginalizing over the foreground model. 

Above about 300 GHz, the DSFGs are the dominant 
source of anisotropy power on all scales. Correlations 
with maps at these higher frequencies, for example maps 
from Planck bands above 217 GHz or Herschel, can be 
used to place tight constraints on the DSFG components, 
at the price of requiring more sophisticated modeling for 
the spectral dependence and shape. Even with such mod- 
eling, the correlations arc no longer fully coherent across 
frequencies so there is a limit to how much of the DSFG 
power can be "clea ned out" of the lowe r frequency maps. 
Following results in Knox et al. (2001 1, which assumes a 
redshift dependent grey- body emissivity density tracing 
the linear matter power spectrum, we assume that we 
could clean out 90% of DSFG clustering power at the 
lower frequencies. As in the previous case of adding in 
Ground data, this again is enough to eliminate nearly all 
of the degradation on cosmological parameters. 



8. CONCLUSIONS 

To make full use of Planck'' s very small statistical er- 
ror on CMB power spectra out to £ ~ 2500, without 
introducing significant bias in the cosmological param- 
eters, we must include contributions from extragalactic 
foregrounds and secondaries in our model of the data. 
Here we have presented a model of these contaminants, 
based on the latest data and modeling developments, 
and demonstrated its ability to remove biases in an 8- 
parameter cosmological model. The foreground model 
has 17 parameters - many more than any extragalac- 
tic foreground model used in analysis of CMB data to 
date. Despite the large number of nuisance parameters, 
marginalizing over all of them only increases statistical 
uncertainties in the cosmological parameters by, at most 
10 to 20%. Almost all of this degradation can be avoided 
by inclusion of ground-based data or higher frequency 
Planck bands. 

Our model includes Poisson components from both ra- 
dio galaxies and DSFGs, a clustering component due to 



DSFGs, contributions to kSZ power from patchy reion- 
ization, as well as after reionization is complete, and tSZ 
power. If kSZ power and tSZ power are at our fiducial 
values (slightly higher than the preferred values given 
current high-resolution ground-based data) then ignor- 
ing them in an analysis of Planck data would produce 
small, almost negligible biases, to cosmological parame- 
ter estimates. On the other hand, ignoring the clustering 
of DSFGs, would lead to a very large bias in cosmological 
parameters. 

To avoid having to marginalize over these 17 param- 
eters every time a new cosmological model is analyzed, 
we broke our procedure up into a two-step process, with 
the first step independent of the model of the primary 
CMB power spectra. The second step is an analysis of 
the CMB power spectra estimated in the first step, with a 
small number of foreground template amplitude param- 
eters to marginalize over. The shapes of these templates, 
and priors on their amplitudes, are also outputs of the 
first step. Only the second step needs to be repeated 
in order to get constraints on the parameters of a new 
model of the primary CMB power spectra. 

Looking toward the near future, the model will def- 
initely evolve, increasing the faithfulness with which it 
represents reality, as we gain more information from the 
CMB-dominated channels in Planck, higher-frequency 
Planck channels, higher-resolution ground-based data 
(SPT, SPTpol and ACTpol) and higher-resolution, 
higher-frequency space-based data (Herschel). One could 
easily use our foreground model to study potential biases 
in extensions of the primary cosmological model, to in- 
clude, for example, departure of the Helium mass fraction 
from its BBN-expected value, or the number of effective 
neutrino species from its standard model value. 
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APPENDIX 



A. CMB LINEAR COMBINATION GENERALIZATION TO OFF-DIAGONAL CORRELATIONS 

The method for constructing a best estimate of the CMB presented in Section [6] assumes only temperature power 
spectra, and a covariance which is diagonal in t. The generalization to include polarization and mode-mode coupling 
induced by sky masking is presented here. The math is, infact, indentical for the two scenarios, so in this appendix 
we'll refer to polarization types with the understanding that we could just as well be talking about different values of 

The added difficulty in dealing with different power spectrum types (for simplicity here just TT and EE) comes from 
the fact that we cannot arbitrarily create linear combinations which sum them. For example, C" — C^^^qyit. ~ ^u)o'ghz 
neither preserves CMB normalization, nor can we be sure it is CMB-free independent of model. To remedy this, we 
make sure that in our construction, any linear combination we consider must have the CMB signal cancel out for all 
but one type. For example, C" = C^'JJghz ~ ^kio'ghz + '-^mcHz ^ valid linear combination. 

We start by considering the covariance matrix for the TT and EE spectra. 



^EE 



(Al) 



By creating the single-type weight matrix (Eq. 



28 1 for each of the diagonal blocks, we can cancel the CMB out of all 
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but two weightings. The new covariance wiU look hke, 



Wf-r 




^rprp . . . 

■ ■ ■ ^EE 



Wtt 
Wee 



( CTEE 

Under a permutation to place the two CMB weightings at the front, the covariance becomes, 




(A2) 



CFrprp . . . 

•■■ CTee 



^cmb '^cross 
'^cross '^diff 



(A3) 



where we've labeled Scmfi as the covariance between TT and EE estimates, Edi// as the covariance of the CMB-free 
differenced spectra, and "Ecross as the cross-correlation between the two. We now would like to do one final reweighting 
in an attempt to zero out the cross-correlation. The reweighting should leave the differenced spectra unchanged, should 
not add TT and EE together, but will add CMB and CMB-free power spectra. Note that this will continue to satisfy 
our earlier condition that all but one CMB type canceling out. The reweighting matrix will look like. 



/ 
W I 



The new covariance must satisfy. 



/ 



^cmb ■^cross 
'Across ^dif f 



I 
W I 



Solving for W yeilds. 



^* — ^fllj ^^cross 



(A4) 

(A5) 
(A6) 



B. PRINCIPAL COMPONENT ANALYSIS FOR POWER SPECTRA 



In this paper we use a p rincipal component analysis (PCA) to reduce the dimensionality of the tSZ as trop hysical 
parameter space (Sec. 2.2) and of the entire foreground contribution to the CMB linear combination (Sec. 6.2). Here 
we present in more detail the procedure used in those sections. 

Given Usim realizations of an n^-length power spectrum, drawn from a statistically significant sample of parameter 
space, we first form the [n^ x nsim] matrix Y . In each column of Y we place the deviation from the mean power 
spectrum for that realization. This matrix is then subject to a singular value decomposition. 



Y = USV'^ 



(Bl) 



where the columns of U contain the orthogonal basis vectors, S' is a diagonal matrix of the singular values and the 
columns of V are the principle component weights. The z-th realization can be written as 



where the singular value-weighted orthogonal basis vectors are 



1 



and the are the weights, 



IfJ, 



(B2) 



(B3) 



(B4) 



Because the singular values are in decreasing order, we can truncate the sum in Eq. |B2| at some small value of ^ 
and still accurately describe each realization. Furthermore, the distribution of weights P(w^) sampled over risim, 
realizations provides a prior on our principal component amplitudes equivalent to the parameter space which was 
sampled to produce the Y matrix. 



